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The nonlocal nature of the fermionic mode spanned by a pair of Majorana bound states in a 
one-dimensional topological superconductor has inspired many proposals aiming at demonstrating 
this property in transport. In particular, transport through the mode from a lead attached to the 
left bound state to a lead attached to the right will result in current cross-correlations. For ideal 
zero modes on a grounded superconductor, the cross-correlations are however completely suppressed 
in favor of purely local Andreev reflection. In order to obtain a non-vanishing cross-correlation, pre¬ 
vious studies have required the presence of an additional global charging energy. Adding nonlocal 
terms in the form of a global charging energy to the Hamiltonian when testing the intrinsic nonlo¬ 
cality of the Majorana modes seems to be conceptually troublesome. Here, we show that a floating 
superconductor allows to observe nonlocal current correlations in the absence of charging energy. We 
show that the non-interacting and the Coulomb-blockade regime have the same peak conductance 
e^//i but different shot-noise power; while the shot noise is sub-Poissonian in the Coulomb-blockade 
regime in the large bias limit, Poissonian shot noise is generically obtained in the non-interacting 
case. 

PACS numbers: 74.45.-l-c, 74.25.F-, 73.23.-b, 74.78.-w 


I. INTRODUCTION 


A pair of Majorana bound states at the ends of one¬ 
dimensional p-wave superconductors hosts a nonlocal 
fermionic mode at an energy close to the middle of the 
energy gap>i The nonlocality of the fermionic mode has 
inspired proposals to observe nonlocal current correla¬ 
tions which arise when electrons are transported through 
the mode from a lead contacting the left bound state to a 
lead contacting the right.— Unfortunately, nonlocal trans¬ 
port is completely suppressed in favor of local Andreev 
reflection for ideal Majorana bound states at zero energy 
and the currents through the left and right lead are un¬ 
correlated with characteristic peak conductances 2e^ jh^ 
Later, it has been pointed out that a coupling t between 
the Majorana end states that is larger then the lead- 
induced level broadening F is required for a recovery of 
nonlocal current correlationsBy increasing the length 
of the superconductor L the coupling constant t decays 
exponentially such that for sufficiently long topological 
superconductors the transport remains local despite the 
nonlocal nature of the fermionic zero mode. 

It was first shown in Ref. @ that to observe nonlocal 
current cross-correlations for ideal Majorana zero modes 
with t = 0, it is important that the superconductor car¬ 
rying the Majorana zero modes is disconnected from the 
ground plate, resulting in a finite charging energy. In 
the Coulomb blockade regime, local Andreev reflection is 
blocked energetically and nonlocal transport with a sys¬ 
tem size independent level broadening F and a peak con¬ 
ductance of /h at zero temperature is obtained. Due to 
the suppression of local Andreev reflection, the analysis 
can be carried out in the usual mean-field setup of su¬ 
perconductivity without violating current conservation.— 
Subsequent studies have considered the behavior of the 
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FIG. 1. (a) Setup consisting of a one-dimensional topo¬ 

logical superconductor (TS) that is floating and hosts ideal 
zero-energy Majorana bound states 71, 72 contacted by two 
normal-conducting, non-interacting leads (N). The supercon¬ 
ductor is subject to a charging energy Ec such that a pre¬ 
ferred number of charges is controllable through a gate volt¬ 
age Vg. We study the system In the case Ec = 0 and compare 
it to known results for large Ec- {b) In absence of charging 
energy, Ec = 0, the conductance and shot noise of the system 
can be understood in terms of an equivalent circuit consisting 
of (noiseless) conductances Gi which are associated with local 
Andreev reflection off the bound state ji and current sources 
generating noise currents SR. 


peak conductance as the charging energy is lowered away 
from the Coulomb blockade regime, finding a smooth in¬ 
crease from the peak conductance e^jh in the Coulomb 
blockade regime to 2e^ jh in the absence of charging 

energy— 2iii 

The system size independence of F in Ref. must 
be contrasted with the transmission of spinless electrons 
through a double barrier, where the level broadening F 
scales with system size as L~^ such that the conduc¬ 
tance quantization at a value jh is lost for any finite 
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temperature in the infinite system limitIn contrast, 
the peak conductance /h survives at finite tempera¬ 
tures in the Major ana case, justifying the terminology of 
calling it ^nonlocal transport'. However, the presence of 
the large global charging energy, which is a strong, non¬ 
local perturbation at the level of the Hamiltonian, makes 
it difficult to unambiguously attribute the nonlocal cur¬ 
rent correlations to the presence of the Majorana bound 
statesJ^ In this sense, a large charging energy can be 
considered ill-suited for probing the Majorana-induced 
nonlocality. 

In this work, we consider a floating superconductor 
where the charging energy Eq vanishes. We will contrast 
the results to the well-studied limit of Coulomb-blockade 
Ec —t oo where the model maps onto a resonant-level 
models We show that consideration of a floating topo¬ 
logical superconductor in a setup as depicted in Fig.jlja) 
allows to observe nonlocal current correlations for ideal 
Majorana bound states at zero energy in the absence of 
charging effects. Interestingly, we find that the peak con¬ 
ductance is given by e^//i both in the non-interacting and 
in the Coulomb blockade limit. The contradiction to the 
results of Refs. M and [ll| can be resolved by noting that 
away from the Coulomb blockade regime, transport into 
the Cooper pair condensate is no longer blocked energet¬ 
ically and current conservation has to be enforced self- 
consistentlyi^ This appears to be technically demanding 
in the interacting setup of finite charging energy. Even 
though the peak conductances does not distinguish the 
two regimes, the transport mechanisms differ strongly. 
In the Coulomb-blockade regime, the Cooper-pair con¬ 
densate is irrelevant and transport proceeds by sequen¬ 
tial tunneling which can be modeled as coherent trans¬ 
port through a double barrierJ^^— In contrast, in the 
non-interacting case, transport is incoherent due to the 
presence of the Cooper-pair condensate which acts as a 
phase-breaking scattererJi We will show that while the 
difference between the two transport regimes is not cap¬ 
tured by the peak conductance, it is reflected in the shot 
noise properties. 

The effective phase-breaking in electron transport due 
to the Cooper-pair condensate allows to understand the 
transport in terms of an equivalent, classical circuit de¬ 
picted in Fig. [TJ&). In this circuit, Andreev reflection 
events at the left and the right terminal are modeled by 
noiseless conductances Gi, G 2 and their corresponding 
noise currents Sli are generated by current sources in par¬ 
allel to the conductances. Straightforward application of 
Kirchhoff rules yields the total conductance G = (/) jV = 
GiG 2 /(Gi-|-G 2 ) since the resistances G~^ simply add up. 
In presence of ideal Majorana bound states, the each of 
the two conductances Gi and G 2 in series peak at a value 
equal to the conductance quantum 2e^/h. Accordingly, 
the total conductance G has a peak vale of half a conduc¬ 
tance quantum e^/h. Similarly, assuming uncorrelated 
noise currents {dli{t)dl 2 {t')) = 0 , one obtains the total 
current noise-power S = Jdt (SI{0)6I{t)) = -I- 

G?4°V(G'i +G 2)2 with = Jdt {SIjiO)SIjit))- here. 


the prefactors in front of the factors originate from 
current division and give the (squared) fraction of the 
noise currents Sit contributing to the total current /. The 
above results correspond directly to our results Eqs. © 
and (IH), yielding an appealingly simple picture of trans¬ 
port through the system. By analyzing the interplay be¬ 
tween the bias dependencies of the conductances and the 
noise currents, we will show that the shot noise in the 
non-interacting limit Ec = 0 generically becomes Pois- 
sonian in the large-bias limit. In contrast, the shot noise 
is sub-Poissonian in the Coulomb-blockade regime. 

The outline of the paper is as follows. We use Sec. Ullto 
fix our notation and give a brief review of results obtained 
in previous studies for the grounded case. We discuss the 
conductance and shot noise obtained in the floating case 
for Ac = 0 in Sec. HH We will compare the results to 
the well-studied case of a superconducting island in the 
Coulomb-blockade regime in Sec. IIVI and finish with a 
discussion of our results. 


II. REVIEW OF THE GROUNDED CASE 

Our system of interest is a one-dimensional p-wave 
spinless superconductor that hosts a pair of Majorana 
bound states 7 ^, i = 1 , 2 , each of them in tunnel con¬ 
tact with a single normal-conducting lead i as depicted 
in Fig. UKa). Before discussing the case of a floating 
superconductor that is disconnected from the ground, 
we start with a review of the grounded case. For the 
grounded case, we describe the system using a mean- 
field Bogoliubov-de-Gennes (BdG) Hamiltonian and cal¬ 
culate the transport properties using a scattering ap¬ 
proach. The entries of the scattering matrix describe 
the transport of particle of type f3 from contact j to con¬ 
tact i as a particle of type a, where i,j S {1,2,5'} can 
be either of the normal-conducting terminals or the su¬ 
perconductor (5) and a, /3 € {e, /ij denote electrons and 
holes. With the transfer probabilities the 

steady-state current at the normal-conducting terminal i 
into the system can be written 



with fj^E = —E) the Fermi distribution of contact j 

at an electro-chemical potential pLj at zero temperatures 
with Q(x) the unit step function. In the following, we 
choose a gauge where the energy E is measured with re¬ 
spect to the chemical potential /ig of the superconductor 
such that fij = eVj with Vj denoting the voltage dif¬ 
ference between the lead j and the superconductor. The 
transport properties are described by the differential con¬ 
ductances Gij = dli/dVj which have the form 

r)T 

( 2 ) 
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at zero temperature; here, T^^{eV) denotes the trans¬ 
mission probability at energy eV. For ideal Majorana 
bound states and low voltages with energy below the gap 
of the superconductor, contributions from crossed An¬ 
dreev reflection {T 12 ), normal transmission (Tfl), and 
quasiparticle transport into the supercondutor (T^g) 
vanishji^ Consequently, the only processes which are pos¬ 
sible are normal- (TJ/) and Andreev-reflection {T^f) such 
that the quasiparticle conservation leads to the sum rule 
Tjj + = 1. As a result, the conductance matrix 

becomes local, = SijGj, with 

o, = (3) 

As shown in App. E explicit calculations of the scat¬ 
tering matrix for ideal Majorana bound states yield 
rj*®(0) = 1, implying resonant Andreev reflection with 
a peak conductance of 2e^//i at zero bias Vj = 0.— 

The nonlocal current correlations can be investigated 
through the zero-frequency shot noise 




dt 


{mm)-{if 


( 4 ) 


defined in terms of the current operator li (t) at terminal 
i. For ideal Majorana bound states in a three-terminal 
NSN setup where the superconductor is grounded, the 
shot noise at each terminal takes the same form as in a 
conventional NS junction with^S 

= X / (5) 

and no current cross-correlation is present)^ 


III. FLOATING CASE WITHOUT 
INTERACTIONS 


In the case of the floating superconductor without in¬ 
teractions Ac = 0, the mean-field BdG formalism which 
neglects the Cooper-pair condensate dynamics cannot be 
applied naively since electron pairs can be converted into 
Cooper pairs, resulting in a net current Ig = —(Ii + 12 ) 
out of the superconducting reservoir. For a floating su¬ 
perconductor in the stationary regime, no net current Ig 
can be drawn from the reservoirThe condition Ig = 0 
can always be fulfilled by varying the sum of the voltage 
Vi -I- V 2 (which corresponds to changing the chemical po¬ 
tential of the superconductors with respect to the leads) 
while keeping the voltage difference V = V 1 —V 2 between 
the two terminals fixed. The superconducting reservoir 
then serves as a voltage probe breaking the phase coher¬ 
ence of the quasiparticles without affecting the current 
balance. 

For a general floating NSN setup, one obtains the 

conductanc o^d^ 


G = 


dV 


_ G11G22 — G12G21 

is=o Gii -l- Gi 2 -|- G21 -f G2 


( 6 ) 


where all the conductances Gy have to be evaluated at 
the voltages Vi corresponding to Ig = 0. For ideal Ma¬ 
jorana bound states, we have G 12 = G 21 = 0 and the 
formula simplifies to 


G = 


G 1 G 2 
Gi + G 2 ’ 


( 7 ) 


which corresponds directly to the conductance that we 
have derived in the introduction for the equivalent circuit 
Fig. [Ii;6). As previously explained, the total resistance 
G~^ is simply determined by the sum of the resistances 
Gi^, G^^ associated with Andreev reflection processes at 
the left and right terminal and the total conductance G 
peaks when the conductances Gi, G 2 attain their max¬ 
imum values at Vi = V 2 = 0. At this point, current 
conservation follows trivially from /i = J 2 = 0 and one 
obtains a peak conductance of jh at zero bias V = 0. 
The peak conductance in the floating case is thus pre¬ 
cisely one half of the peak conductance expected for a 
grounded superconductor. 

For the calculation of the noise in the floating case, we 
follow the approach of Ref. which assumes that the cur¬ 
rent fluctuations can be modeled as Langevin forces 51 j 
that drive the fluctuations —e5V in the chemical poten¬ 
tial of the superconductor which in our gauge corresponds 
to modifying the voltages according to dUi /2 = 5V . Al¬ 
though Ref. H presents the results in linear response 
in the voltages Vj where they can be rigorously justi¬ 
fied through fluctuation-dissipation relations,— the treat¬ 
ment formally only requires linear response in the fluc¬ 
tuations 5V of the voltages Vj. This can be seen by 
expanding the zero-frequency current operator Ii of lead 
i in the form 


ii — Ii 5- 'y ( Gij SV + Sii, (8) 

jits 

where the first term is just the average current Eq. (HD 
and 5Ii are the Langevin forces defined by {5Ii{t)) = 0, 
Jdt {Sii{0)6Ij{t)) = with the shot noise S^j'^ of the 
grounded case. In equation ([5]), the expressions Ii, Gij 

and s{^j'^ all have to be evaluated at the voltages corre¬ 
sponding to Ig = 0. For long measurement times which 
corresponds to low frequencies, current conservation has 
to hold not only on average but as an operator identity 

Using the current conservation, the result of Ref. |9| for 
the shot noise in a floating two-terminal setup is recov¬ 
ered. The shot noise is maximally cross-correlated with 
S'!! = S 22 = —S 21 = —S 12 = S where 

ri2 0 ( 0 ) I f^2 q(0) 

U2J1 -|-(-Ti02 /n^ 

* - (G, + G,y ■ 

Note that the expression of the shot noise (O can directly 
be understood in terms of the equivalent circuit Fig.[T](6) 
that was discussed in the introduction. However, in order 
to determine the concrete bias dependence of the noise 
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FIG. 2 . Relation between the voltage drop Vhiin at the termi¬ 
nal with the smaller lead coupling Finin and the voltage drop 
Rmax at the terminal with the larger lead coupling Fmax for a 
floating superconductor. 

(0), we need to solve the self-consistency equation Jg = 0 
determining the voltages Vi given the bias V. We deter¬ 
mine this relation using a microscopic model of two ideal 
Majorana bound states 71, 72 that are each connected to 
a lead. Because there is no coupling between the ideal 
Majorana modes, the transport is local and we can use 
the results of App.|A]in the expressions Eqs. o, © and 
(IH) for Ij, Gj, and We obtain 

Gj = ^(1 + )i Ij = arctan(wj), (10) 

(arctan \vj \ - \vj \/(1 -h i;|)), (11) 

with Vj = eVj/ 2 Tj where we have introduced the cou¬ 
pling strength Fj of the Majorana fermion 7^ to the 
nearby lead. 

For the following, a possible asymmetry between the 
lead couplings F^- will be important. Thus, we will label 
the leads by the indices j S {min, max} such that Fmin < 
Fmax. Note, however, that for the voltage drop t4iin, 
the label ‘min’ merely means that the voltage drops at 
the terminal with the smaller lead coupling Fmin which 
implies that Vmin > Knax as we will show below. 

Solving the constraint Ig = 0 valid for a floating su¬ 
perconducting in terms of Umax yields 

Umax — tanj^Fmin arCtan(Umin)/rmax] ■ (1^) 

The physics behind the transport can be entirely under¬ 
stood through the voltage relation (TT^ , which is plotted 
in Fig. [^for different ratios Fmin/Tmax- We find that for 
small applied bias voltages F <C F_ = 2Fmin/e the volt¬ 
age drop is symmetric with | 14 iax| = IKninl = 5^- This 
results corresponds to the linearized regime in Fig. [ 5 ] 
Note that for the special case of symmetric couplings this 
relation holds exactly for arbitrary voltages. At larger 
voltage bias, the voltage Fnax saturates at a value 

Fat = 2Fitan(7rFmin/2Fmax)- 



F/Fat 

FIG. 3 . Fano factor F = S/el in the floating sys¬ 
tem as a function of the voltage F/Fat, where eFaat = 
2rmax tan(rniinn'/2rmax) is the largest voltage Fmax that can 
drop at the terminal with the larger lead coupling Fmax. For 
voltages F S> Fsat, nearly all the voltage drops at the terminal 
with the smaller lead coupling Fmin and both the shot noise 
S/e and the current I tend to the maximal current 2 'KeTniin/h 
that can be driven through the system. Consequently, for 
voltages F S> Fat , the Fano factor reaches one. 


The origin of the saturation is the effect that a single An¬ 
dreev level can carry at most a current 27reF/ h. The total 
current through the device is thus limited by 27reFmin/^- 
In order for Is = 0 to hold, the current through the junc¬ 
tion with the coupling Fmax must not exceed this value. 
Using (I10|) . with /max = 27reFmin/^, we obtain that the 
voltage Umax has to saturate at the value Fat- If the 
bias voltage exceeds this saturation value almost all the 
voltage drop will be across the junction with the smaller 
coupling, i.e., we have F = Fmin for F Fat- We high¬ 
light that the existence of this regime is a consequence 
of the resonant structure of Andreev reflection at the NS 
interfaces. 

From the behavior of the voltages, all the other trans¬ 
port characteristics follow. In the regime of low bias 
F <g; F_, the relation have Fnin = — Fnax leads to a non¬ 
trivial relation between the conductances Gmin/Gmax = 
1 - (Tmin + rmax)e^l^Vl6r^i„F^,^,, and the behavior of 
the shot noise in the floating case is determined both by 
the bias dependence of the conductances and the shot 
noises Smax- For F <C F_, we obtain a Fano factor 
F = S/el = {1 + FF„/F^3^^)FV24F2 < 1 which in¬ 
dicates sub-Poissonian noise. For F 0, we obtain the 
result F — 0 consistent with the fact that the local trans¬ 
port is carried by a resonant level with unit transmission 
probability implying a vanishing Fano factor due to the 
Pauli exclusion principleJ^ 

Since for F » Fat, the voltage drop Fnin at the ter¬ 
minal with the smaller coupling Fmin is much larger than 
the voltage drop Fnax at the other terminal, the conduc¬ 
tances in that regime obey Gmin “C Gmax- Consequently, 
we hnd S = which reflects the fact that accord¬ 
ing to the circuit Fig. [T}/), it is the noise current gen- 
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erated at the terminal associated with Finin which flows 
through the larger conductance and thus dominates the 
shot noise for V ^ t4at- Using the explicit form of the 
noise dm and the current (nni), one finds that the Fano 
factor F = Sjel tends to one for V ^ Usat- This behav¬ 
ior is illustrated in Fig. [3] for different ratios Fmin/Fmax 
close to one. 


IV. COMPARISON TO THE 
COULOMB-BLOCKADE CASE 

The behavior of the shot noise power ([9|) in the non¬ 
interacting case Ec = 0 must be contrasted with the shot 
noise power obtained for the Coulomb-blockade regime 
Ec —)• oo of Ref. d. In the case where only two charge 
states are important, transport proceeds by sequential 
tunneling and can be mapped onto transport though 
a double barrier where tunneling through the barriers 
occurs at a rate Fniin/fi, Fmax/fi-— Current and shot 
noise in that case can be obtained using the classical 
circuit Fig. mb) and its corresponding shot noise ex¬ 
pression if the conductances Gmim Gmax model the 
transmission through a single barrier.— In this case, the 
conductances Gmin, Gmax are simply proportional to the 
couplings Fmin, Tmax which are energy-independent for 
wide leads, showing that the voltage drops Umin, Umax 
at the left and the right barrier must be distributed ac¬ 
cording to Umin/Umax = Tmax/Tmin- Since neither of the 
voltage drops saturates, both the Fano factors associ¬ 
ated with transmission through the left and right barrier 
s[^'^/el, /el will tend to one in the large bias regime 
eV ^ Fmin, Fmax and one obtains the limiting expression 
Fdh = (FU„-bF^ 3 ^x)/(Fmin+Fmax)^ < 1 for the Fano fac¬ 
tor of the double barrier in the large bias limit This 
shows that the Fano factor is suppressed as compared to 
the single barrier case with a maximal suppression of 1/2 
obtained for symmetric coupling Fmin = Fmax- 

In the following, we want to give an intuitive picture for 
the reduction of the Fano factor in the Coulomb-blockade 
regime with respect to the case of a noninteracting float¬ 
ing superconductor with Ec = 0. In the former case, 
the Fano factor is given by the generic result F^h for 
an asymmetric double barrier which is smaller than one. 
In the latter case, for large bias, the voltage generically 
drops solely over the junction with coupling Fmin and 
the shot noise as well as the current can be evaluated in 
an NS-setup involving only this junction. The transport 
through the NS junction involves Andreev reflections. It 
is a well-known fact that the Andreev reflection processes 
in an NS-junction can also be mapped onto the problem 
of transmission of electrons through a symmetric dou¬ 
ble barrier with the concreted identification explained in 
Fig. m It is however important to realize that the ele¬ 
mentary processes in the transport involve Cooper pairs 
with charge 2e such that the Fano factor is a twice large 
than the one for the symmetric double barrier. Given the 
expression F^h above, we thus obtain the result that the 



FIG. 4. This plot shows the equivalence of the tunneling 
through a symmetric double barrier to the Andreev reflec¬ 
tion off a NS interface, (a) The tunneling amplitudes for 
a double barrier t = Am is a sum of the ampli¬ 
tudes Am = rL)'^tR for m round trips between 

the barriers with the amplitudes Tr/r and Ir/r of reflec¬ 
tion and transmission at the left/right barrier and the phase 
accumulated during one traversal of the barrier region. 
(6) The tunneling amplitude for Andreev reflection at a 
NS interface is similarly given by the sum of the amplitude 
Am = tee'^{rhe'^rhe'^)'^th for m round trips with the ampli¬ 
tudes and for reflection and transmission of an elec¬ 
tron/hole at the tunneling barrier and the amplitudes e*^, 
for Andreev reflection as a hole or electron at the su¬ 
perconductor. For energies much smaller than the gap, the 
particle-hole symmetry dictates that X -|- X = Thus, 

we obtain the resulb^^ that Andreev reflection can be un¬ 
derstood as transmission through a symmetric double barrier 
with Tr = = Tr = if we identify (p = —tx/2 and 

te/h = th/R, = Tr/r. 


Fano factor of the noninteracting setup tends to one for 
large bias. 


V. CONCLUSIONS 


In this work, we have studied two-terminal transport 
between two normal-conducting metallic leads contacting 
the Majorana bound states of a floating topological su¬ 
perconductor in absence of charging energy. The nonlo¬ 
cal current correlations reflecting Majorana-assisted non¬ 
local electron transport that were previously shown for 
the Coulomb-blockade regime^ are also obtained in the 
non-interacting case. The presence of the phase-breaking 
Cooper-pair condensate allows modeling the transport 
processes as an incoherent combination of Andreev re¬ 
flection events at the left and the right NS junction. Us¬ 
ing this model, we have shown that the peak conduc¬ 
tance e^//i previously obtained for the Coulomb-blockade 
regime also follows in the non-interacting case as a sim¬ 
ple consequence of current division. We have shown that 
the difference between the non-interacting case and the 
Coulomb-blockade regime is reflected in the shot noise 
power, which generically becomes Poissonian in the high 
bias limit of the non-interacting case, whereas it stays 
sub-Poissonian in the Coulomb-blockade regime. We 




















6 


thus find that the two transport regimes can be distin¬ 
guished by their shot noise properties. Our findings for 
the peak conductance make it tempting to speculate that 
the peak conductance will actually remain at value e^//i 
for any finite value of the charging energy which however 
remains a problem for future studies. 


puted through the Weidenmiiller formula for the scatter¬ 
ing matrix at energy 


s{u:) = 1 — 2mW^ 


1 

ul2 + iT:WW^ 


W, 


(Al) 
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where the vector W = (w,w*) describes the coupling 
of the Majorana bound state to the electron and hole 
degrees in the lead which is assumed to have unit density 
of states. The scattering matrix can be written as 


s = 



fl + A A \ 

A 1 + aJ ’ 


(A2) 


Appendix A: S-matrix for ideal Majorana bound 
states 

The scattering matrix for an ideal Majorana bound 
state at a NS-interface can be straight-forwardly com¬ 


with A = (iuj/2r — 1)“^ where we have introduced the 
Fermi’s golden rule coupling strength T = 27r|wp to the 
lead. From this one obtains the Andreev reflection prob¬ 
ability T'*® = (1+ uj^/4r^)-\ 
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